Magnetic Fields and Ionized Gas in the inner Galaxy: An Outer 
Scale for Turbulence and the Possible Role of H II Regions 
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ABSTRACT 

We present an analysis of rotation measure (RM) fluctuations from the Test Region of the 
Southern Galactic Plane Survey (SGPS), along with emission measure (EM) fluctuations in the 
same field taken from the Southern H-Alpha Sky Survey Atlas. The structure function of RM 
fluctuations shows a relatively steep slope at small scales (1-5 arcmin), a break in slope to a 
flatter structure function at intermediate scales (5 - 60 arcmin), and a systematic variation of 
the strength of fluctuations as a function of position angle on the sky at the largest scales (60 - 
200 arcmin). The structure function of EM fluctuations shows similar behavior, although the 
lower resolution of the data prevents detection of a possible break in the spectrum. We interpret 
the anisotropy in RM/EM structure on large scales as resulting from a large-scale gradient in 
electron density (and possibly magnetic field) across the region. The break in the slope of the 
RM structure function at scales of ~ 5 arcmin can be explained by contributions from two 
spatially distinct magneto-ionized screens, most likely in the Local and Carina spiral arms. The 
observed structure function then implies that the outer scale of RM fluctuations in these screens 
is ~ 2 pc. Such behavior is in striking contrast to the expectation that interstellar turbulence 
forms an unbroken spectrum from kpc down to AU scales. We conclude that we have identifled 
an additional source of enhanced turbulence, injected on scales of a few pc, possibly seen only 
in the Galactic plane. The most likely source of such turbulence is individual Hn regions from 
relatively low-mass stars, whose characteristic scale size is similar to the outer scale of turbulence 
inferred here. These sources may be the dominant source of density and velocity fluctuations in 
warm ionized gas in the Galactic plane. 



Subject headings: 
radio continuum; 



Galaxy: structure 
ISM — turbulence 



H II regions — ISM: structure — techniques: polarimetric 



^Harvard-Smithsonian Center for Astro- 
physics, 60 Garden Street, Cambridge MA 01238; 
mhaverkorn@cfa.harvard.edu, bgaensler@cfa.harvard.edu 

^Austraha Telescope National Facility, CSIRO, P.O. 
Box 76, Epping, NSW 1710, Australia; naomi.mcclure- 
grifHths@csiro.au 

3 Bolton Fellow 

* Department of Astronomy, University of Minnesota, 
116 Church Street, SE, Minneapolis, MN 55455; 
john@astro.umn.edu 

^School of Physics, University of Sydney, NSW 2006, 
Australia; agreen@physics.usyd.edu.au 



1. INTRODUCTION 

The evidence for the presence of turbulence in 
the interstellar medium (ISM) is overwhelming. 
Although some studies of the characteristics of this 
turbulence (viz. the inner and outer scales, shape 
of the spectrum and power law spectral index) 
indicate the presence of standard incompressible 
Kolmogorov (1941) turbulence, many observations 
indicate different kinds of turbulence, drivers and 
environments. Armstrong, Rickett & Spangler 
(1995) compiled observations of (among others) 
interstellar scattering of pulsars and extragalactic 
sources, dispersion measures of pulsars and rota- 
tion measures (RM) of extragalactic sources in one 
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power spectrum. The result was the so-called "big 
power law in the sky", a Kolmogorov-like power 
spectrum of electron density fluctuations over 12 
orders of magnitude, from a fraction of an AU to 
kiloparsecs. Although the big power law in the 
sky is very well-determined on small scales (up to 
^ 0.001 pc), the extension towards larger scales 
is based only on two types of measurements: the 
RM of extragalactic sources, which gives an upper 
limit to the amount of structure in the electron 
density because the contribution of the magnetic 
field is unknown; and from velocity measurements 
of H I, making assumptions about the connection 
between neutral and ionized material. 

The dominant source for the large-scale energy 
input in the Galaxy is generally assumed to be su- 
pernova explosions (Spitzer 1978; VoUmer & Beck- 
ert 2002; Korpi et al. 1999). Other possible sources 
are superbubbles, massive H ii regions and massive 
stellar winds (Norman & Ferrara 1996), Galac- 
tic fountains, chimneys, or gravitational scattering 
by transient spiral waves (sec Sellwood & Balbus 
(1999) and references therein), gravitational insta- 
bilities in a shearing disk (Elmegrecn, Elmegreen 
& Leitner 2003), or magneto-rotational instabili- 
ties (Mac Low & Klessen 2004). The big power law 
in the sky certainly suggests the classical scenario 
of turbulent energy input on large scales, which 
cascades down to smaller scales until the energy 
is dissipated on the smallest scale (Kolmogorov 
1941). 

However, there are indications of other types of 
turbulence in the ISM as well. Minter & Spangler 
(1996) found a break in the structure function of 
RMs of extragalactic sources and emission mea- 
sures (EM) of the warm ionized gas on scales of a 
few parsec, which they interpreted as a transition 
from three-dimensional to two-dimensional turbu- 
lence as one moves from small to large scales. Ana- 
lytic theory (Goldreich & Sridhar 1995, 1997) and 
simulations (e.g. Cho, Lazarian & Vishniac 2002; 
Maron & Goldreich 2001) show that magnetic 
fields can lead to anisotropic turbulence, which 
is predicted to exhibit the same power law spec- 
tral index as Kolmogorov turbulence. Further- 
more, observations of RMs of extragalactic sources 
(Simonetti & Cordes 1986; Spangler & Reynolds 
1990; Clegg et al. 1992) show higher amplitudes 
of structure in RM in the Galactic plane than out 
of the plane, suggesting the existence of an addi- 



tional source of structure on parsec scales. Re- 
sults from interstellar scattering show a similar 
enhancement of structure in the inner Galaxy, but 
on much smaller scales (Rao & Ananthakrisnan 
1984; Dennison et al. 1984; Anantharamaiah & 
Narayan 1988). 

Structure in the ionized ISM can be studied by 
way of structure functions (SFs). Earlier determi- 
nations of SFs of RM using extragalactic sources 
or the diffuse synchrotron background yielded re- 
sults that varied widely with scale and area in 
the sky. Flat structure functions (indicating no 
structure on the probed scales) found by Simon- 
etti, Cordes & Spangler (1984) on scales larger 
than 4° towards the Galactic pole were interpreted 
as showing structure intrinsic to the extragalac- 
tic sources with a negligible Galactic contribution. 
However, this cannot explain the flat SFs from 
low-latitude extragalactic sources found by Simon- 
etti et al. (1984) on scales > 4° and by Clegg 
et al. (1992) on scales > 1°. Furthermore, Sun 
& Han (2004) have recently found shallow SFs 
of RM in the Galactic plane, and a fiat SF at 
the North Galactic Pole, also from extragalactic 
sources. Haverkorn, Katgert & de Bruyn (2003a) 
studied the structure function of RM from diffuse 
radio emission, and found very shallow slopes for 
two fields at intermediate latitudes. 

In this paper we study the turbulent struc- 
ture in ionized gas in the inner Galactic plane, by 
means of SFs of RM from the diffuse synchrotron 
background, and of SFs of EM from Ha emission, 
both in a region in the inner Galactic plane. 

In Section 2 we present our radio polarimet- 

ric observations of the Galactic synchrotron back- 
ground in the Galactic plane. Section 3 discusses 
the computation of the structure function, while 
in Section 4 we interpret the structure function as 
arising from Faraday screens in two spiral arms 
along the line of sight, both exhibiting turbulent 
structure. In Section 5, we discuss evidence for en- 
hanced density fiuctuations in the Galactic plane, 
and propose that the structure is dominated by 
discrete H ii regions. 

2. THE SOUTHERN GALACTIC PLANE 
SURVEY TEST REGION 

The Southern Galactic Plane Survey (SGPS) is 
a radio survey in Hi and linearly polarized con- 
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Fig. 1. — Rotation measure in the SGPS Test Region in squares, overlaid on polarized intensity in grey scale 
given in Jy/bcain. Open squares denote negative RMs and filled squares positive ones. The length of a square 
is proportional to the magnitude of the RM for |RM| < 100 rad m~^, and constant for |RM| > 100 rad m~^. 
RM values have been given only if S/N > 5 and reduced Xr < 2- The rectangular box drawn in the figure 
shows the area over which structure functions are computed. 
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Fig. 2. — Ha emission in the SGPS Test Region, from the Southern H-Alpha Sky Survey Atlas (Gaustad et 
al. 2001). The scale is in deciRayleighs (1 Rayleigh = 10^/47r photons cm~^ s~^ sr~^) and the drawn box is 
the same as in Figure 1. A gradient in emission is visible from the bottom left corner to the top right corner 
of the box. 
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timmm at 1.4 GHz performed with the Austraha 
Telescope Compact Array (ATCA) and Parkes 
64m dish (McClure-Griffiths et al. 2001; Gaensler 
et al. 2001). The first phase of the survey ex- 
tends over 253° < I < 358° and |6| < 1.5°, has 
a spatial resolution of 1' and a velocity resolu- 
tion of 1 km s~^. Before embarking on the com- 
plete survey, a Test Region was observed, span- 
ning 325.5° < I < 332.5° and -0.5° < b < 3.5° 
(Gaensler et al. 2001). The continuum data for 
the SGPS Test Region consists of nine frequency 
bands each of 8 MHz in width, in the range 
1336 MHz to 1424 MHz. From the hnear po- 
larization position angle (j) at these nine frequen- 
cies, RMs were derived as (/> oc RMA^, where 
RM = 0.81/ n^B\\ds with the thermal elec- 
tron density in cm~"^, the component of the 
magnetic field parallel to the line of sight in /xG, 
and ds the path length in parsecs. Gaensler et al. 
(2001) give a detailed analysis of the structure in 
linear polarization and RM in the Test Region, fo- 
cused on morphology and individual sources. RM 
fluctuations were apparent in their Figure 7, but 
not analyzed. We use their data, but use a slightly 
different algorithm to derive RMs. 

Gaensler et al. use a standard algorithm within 
the data reduction package MIRIAD^ (Sault & 
Killeen 2003) to compute RMs, which they de- 
scribe in detail. The MIRIAD algorithm solves 
the problem of wk radians ambiguity in the polar- 
ization angle by flrst computing an RM value from 
the first two frequencies. Then, the polarization 
angles of the other frequencies are rotated titt ra- 
dians so that they lie closest to the angle predicted 
from the first two frequencies, and the RM is re- 
computed using these angles. Therefore, a large 
error in (one of) the two first frequencies can yield 
the wrong determination of RM. In this way, a 
small part (~ 6%) of the RM values Gaensler et 
al. computed was probably incorrect. Note that 
all these erratically determined RMs had high re- 
duced of the relation cj) = RMA^, which caused 
them to be disregarded in their analysis. There- 
fore the use of this algorithm did not lead to wrong 
quoted values of RM, but rather to a small number 
of unjustly disregarded RMs. 

Here, we incorporate the ambiguity of polar- 
ization angle (j) by always taking the least pos- 

^ http ; / / www. atnf . csiro.au/ computing/software /miriad 



sible angle difference between in adjacent fre- 
quency bands, i.e. — 4>j\ < tt rad where i 
and j are adjacent frequencies (Haverkorn, Kat- 
gert & de Bruyn 2003b; Brown, Taylor & Jackel 
2003). The frequencies are so closely spaced that 
this method renders the correct RM for all RM < 
1900 rad m~^. Since in our field RMs are not 
higher than a few hundred rad m~^ and mostly 
much lower, this is the case for all derived RMs. 

Figure 1 shows the RM in the Test Region, su- 
perimposed on the continuum polarized intensity 
in grey scale. The resolution is about 1' and the 
data are about three times oversampled. Open 
(filled) squares represent negative (positive) RM, 
and the size of a square is proportional to the mag- 
nitude of RM for \RM\ < 100 rad m'^ and con- 
stant for \RM\ > 100 rad m~^. Only every third 
beam is plotted for clarity. Furthermore, only "re- 
liably determined" RMs are shown, i.e. RMs for 
which the signal to noise ratio S/N > 5 and the 
reduced of the linear (^(A^) fit was Xr < 2. The 
dependence of (f) on A^ can be non-linear due to 
depolarization effects (Burn 1966; Sokoloff et al. 
1998). Mild non-linearity is allowed by including 
reduced Xr up to 2. About 41% of the data have 
low enough reduced Xr^ of which 37% has a high 
enough S/N. So about 15% of the data has reliably 
determined RMs. 

As these data are taken solely with the ATCA 
interferometer, they lack large-scale information 
(> 30') in the observed Stokes Q and U maps. If a 
large component of Q and U is missing, the com- 
puted values of polarization angle will be incor- 
rect, and therefore the RM will be in error as well. 
However, the polarization angle depends very non- 
linearly on Q and U, so that a large missing com- 
ponent in Q and/or U would destroy the linear 
relation between <+> and A^ . The fact that so many 
RM's with low reduced are observed suggests 
that any missing large-scale components should be 
small. This needs to be confirmed with data from 
the single dish Parkes telescope, which will be the 
subject of future studies. 

We note that the inability to detect large-scale 
structure in Stokes Q and U images due to the lack 
of single-dish data does not imply an inability to 
detect large-scale components in RM. As long as 
small-scale structure in RM exists superimposed 
on any large-scale RM, the complete Q and U sig- 
nals will be detected, which correspondingly allow 
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us to retrieve the total RM, not just the small-scale 
component. Therefore, fluctuations in RM can be 
probed at much larger scales than the maximum 
spatial scale probed by the interferometer. 

As RM only provides information on the prod- 
uct nei3|| integrated over the path length, inde- 
pendent information on the path length integrated 
distribution of thermal electrons is extremely use- 
ful for estimation of the structure in the mag- 
netic field. Therefore, we also used emission mea- 
sure EM = J n^ds data from the Southern El- 
Alpha Sky Survey Atlas (SHASSA, Gaustad et 
al. 2001), which provides comparable resolution 
Ha data of the Southern sky, see Figure 2. Cau- 
tion is needed in comparing RM and EM, as the 
path lengths may not be the same, and the Ha 
data is not corrected for dust extinction. The 
SHASSA provides data which is median filtered 
over four 0.8' pixels to remove star residuals, and 
has a sensitivity of about 0.5 Rayleigh (1 R = 
10^/47r photons cm~^ sr~^, corresponding to 
EM w 2 cm-6 pc at T = 8000 K). 

3. DETERMINATION OF STRUCTURE 
FUNCTIONS 

Ideally, the scales of variability in a signal such 

as RM are determined by the power spectrum. 
However, as the RM data contain blank spaces in 
the regular grid, we use structure functions (SFs) 
instead (Simonetti et al. 1984; Minter & Spangler 
1996). We define the second order structure func- 
tion SF of RM and EM as a function of distance 
lag r as 

SFrmW = ((RM(x)-RM(x + r))2), (1) 
SFemW = ((EM(x)-EM(x + r))2), (2) 

where ()x denotes the averaging over all positions 
X in the field. 

Structure functions of RM and EM were de- 
termined in the area of the SGPS Test Region for 
which the most data were available, i.e. within the 
box drawn in Figs. 1 and 2. The noise was taken 
into account by subtracting a SF of the error 6 
in the RM data from SFrm, as explained in the 
Appendix. We computed one-dimensional SFs for 
different position angles a (positive b through pos- 
itive I) of the distance lag vector r, with an interval 
of Aa = 10°. 



The left hand plots in Figure 3 give the SFs 
of RM (top) and EM (bottom) for different po- 
sition angles 0° < a < 180° superimposed for 
regular sampling in r. Errors in the SF are typ- 
ically log(SFRM) ^ 0.05. The two SFs agree re- 
markably well: both are very shallow; the slope of 
SFrm is ^ 0.2 and the slope of SFem ^ 0.25 for 
4' < r < 5°. The slope of SFrm shows a break, be- 
low which the SF steepens to ~ 0.5. Furthermore, 
at larger scales an anisotropy occurs with position 
angle in both SFs. A linear fit to the slope at large 
scales was performed for each position angle over 
the range given by the horizontal lines in the left 
hand plots. The slopes of these linear fits are given 
in the right hand plots of Figure 3. The anisotropy 
on large scales in the SFs of both RM and EM is 
more or less sinusoidal, with a maximum for RM 
and EM at the same position angle a « —50°. 

The sudden change in slope on scales smaller 
than the break is probably not due to resolution 
effects in the RM data, as the resolution of the 
RM data is FWHM = 86" (log(r) « -1.6), which 
is well below the break. In the EM data, a me- 
dian filtering over 4 pixels means that at scales 
log(r) > —1.17, beams are independent, so that 
the position of the break in SFrm falls below the 
resolution of the Ha data. The unsmoothed data 
(with a resolution of 0.8') are too contaminated 
with point sources to determine if a break in the 
spectrum is present in those data. 

The anisotropy in the SF slope at large scales 
is not an artifact of the gridding or the coordinate 
system, as its maximum does not coincide with 
one of the axes. Furthermore, we constructed one- 
dimensional SF's of the same data in (RA. dec) 
coordinates. The slopes of these SF's show max- 
ima at the same intrinsic position angle as the 
data in Galactic coordinates, confirming that the 
anisotropy is not an artifact of the gridding of the 
data. The same analysis has also been applied to 
maps of the Stokes Q and U parameters, and to a 
map of randomly scrambled RM values in which 
the observed RMs are redistributed at random po- 
sitions in the field. The SF's of none of these maps 
showed an anisotropy, again indicating that the 
change in slope of SF with direction is physical 
and not an artifact. 
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Fig. 3. — Top left: Structure functions of RM, SFrm, in the SGPS Test Region for different position angles 

0° < a < 180° (positive b through positive I) superimposed. The vertical dotted line shows the position of 
the break in slope, and the horizontal line gives the range over which the linear fit was performed to yield 
the slope given in the right hand plot. Top right: dependence of the slope of SFrm on position angle a. 
Bottom: The same for SFem- 
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4. INTERPRETATION OF THE STRUC- 
TURE FUNCTIONS 

4.1. The Anisotropy of the SF Slope at 

Large Scales 

Two very difFcreiit characteristics of a medium 
can produce an anisotropy of the slopes of SFs: 
anisotropic turbulence and a large-scale gradient 
in the medium. 

If a strong enough large-scale magnetic field 
is present in the turbulent ionized medium, the 
turbulence will be anisotropic (e.g. Goldreich & 
Sridhar 1995, 1997; Cho & Vishniac 2000). Ob- 
serving anisotropic turbulence through RM data 
is believed to be difficult or even impossible be- 
cause the anisotropy is directed with respect to 
the local magnetic field (Cho et al. 2002), so that 
structure in the magnetic field along the line of 
sight is likely to destroy any observational effects 
of the anisotropic turbulence. If anisotropic tur- 
bulence is present in the SFs, it should be seen on 
all scales instead of only large ones. Furthermore, 
anisotropic turbulence would yield a set of SFs as a 
function of position angle which have the same am- 
plitude of fluctuations but on different scales. This 
is not the case in our observations, where slopes 
vary from ~ 1 to below zero with position angle. 
Therefore, we do not believe that the anisotropy 
in the slope of the SF is caused by anisotropic tur- 
bulence. 

If the anisotropy in the slope is caused by a 
gradient across the field, one would expect a steep 
SF slope in the direction of the gradient, and a 
flat SF slope in the perpendicular direction. Noise 
and the superposition of a regular structure like a 
sinusoid can decrease the slope to below zero. This 
variation in slope to values below zero is exactly 
what we observe in both SFrm and SFem- In 
fact, in the Ha data the gradient can be directly 
seen in the maps (Figure 2). An Ha map of a 
larger region shows a complex of H ii regions and 
supernova remnants extending over more than 15 
degrees in Galactic longitude at higher longitude 
and lower latitude from the SGPS Test Region, 
and the EM from that complex decreasing towards 
and across the Test Region. This coincides with 
the direction of the gradient computed from the 
SFs. 

A maximum amount of fluctuations in RM at 



the position angle —50° docs not determine if the 
gradient is directed towards —50° or 130°. To find 
the direction of the RM gradient, we looked at 
the correlation between the magnitude of EM and 
RM. Figure 4 shows the average EM in bins of 
|RM| with width ARM = 10 rad m'^. The mag- 
nitudes of EM and RM arc clearly correlated, in- 
dicating that the gradient in |RM| follows the gra- 
dient in EM. Therefore |RM|, like EM, decreases 
towards a = —50°. This is to be expected if gradi- 
ents in electron density and/or path length domi- 
nate, and the magnetic field either shows the same 
gradient or does not change significantly over the 
field. 

To determine what kind of gradient could pro- 
duce an anisotropy in the slope of SF as observed, 
we modeled the anisotropy in SFrm using three 
simple functions: a linear gradient, a broad sinu- 
soid and a Gaussian decrease, all with the gradient 
oriented at a = —50°. Figure 5 shows the results 
of the modeling with each of these three functions, 
compared to the data given by the diamonds. The 
Hnear gradient is given by the dotted line, the si- 
nusoid by the dashed line, and the Gaussian is de- 
noted by the dash-dotted curve. Noise is included 
in the models at RMndse = 5 rad m~^, similar to 
the noise level in the data. These simplified mod- 
els show an anisotropy in the slope of comparable 
amplitude and minimum to the data, although we 
cannot distinguish between them. More elaborate 
modeling of the anisotropy in the SF slope is be- 
yond the scope of this paper. 

The three models fit best for a total gradient 
in RM across the field of ARM « 3 rad m'^, 
while the Ha data show a change in EM of about 
10 cm-f' pc (- 5 Raylcigh at T s=: 8000 K). These 
changes can be duo to gradients in electron den- 
sity or the path length, and most likely a combina- 
tion of these, possibly accompanied by a gradient 
in magnetic field. Therefore the measured ARM 
and AEM give too few constraints to estimate the 
gradients separately in magnetic field strength and 
electron density. 

4.2. The Shallow Slope and Break 

A structure function of RM produced by a slab 
of Kolmogorov turbulence with an inner scale r, 
and an outer scale Tq is shown as the solid line 
in the upper panel of Figure 6 (e.g. Stinebring 
& Condon 1990). The SF is constant on small 
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Fig. 4. — Average EM in bins of absolute RM with 

width AjRMl = 10 rad m~^. The errors are the 
errors in the mean within a bin. 
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Fig. 5. — Slope of SFrm as a function of posi- 
tion angle a (diamonds), and models of a linear 
gradient (dotted), a broad sinusoid (dashed), and 
a Gaussian (dash-dotted). A noise component of 
5 rad m~^ has been added to the model RMs, 
comparable to the noise in the data. 
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Fig. 6. — Example of the superposition of two 
Kolmogorov-like structure functions. Top: two 
structure functions of turbulent Faraday screens 
with an outer scale and an inner scale r.,;, at 
distances "near" (dashed) and "far" (solid). Bot- 
tom: the addition of the two structure functions 
in the top panel shows a break in the slope, with 
a shallower slope at large scales. 
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scales due to dissipation at the inner scale, and 
then starts to rise at rj. If the spatial fluctu- 
ation in RM is a statistically homogeneous and 
isotropic random variable and the spectrum of RM 
is a power law with spectral index s, the slope of 
the SF b (SFrm oc r** for ri <r < ro) is related to 
s as (Simonetti et al. 1984): 

CTT /\ b - , ( s-2 2<s<4 
SFRMi^-jocr where 6 ^ 

(3) 

Therefore a Kolmogorov spectrum (s = 11/3) 

should exhibit a structure function slope b = 5/3. 
Above the outer scale of structure Tq, the structure 
function saturates to the constant value of 2cr^]y[, 
where (t^^;[ is the variance of the RM distribution. 

However, the observed SFs do not show this 
behavior. In the observed SFrm, a break is visi- 
ble, and the slope at scales larger than the break 
is very shallow. This behavior can be most eas- 
ily explained by the superposition of two Fara- 
day screens at different distances, both with a 
Kolmogorov-like SF and with certain physical in- 
ner and outer scales. Even if the physical scales of 
the structure in these screens are equal, the RM 
will exhibit structure on different angular scales, 
as shown in Figure 6. The upper panel shows two 
Kolmogorov-like SFs with the same amplitude at 
different angular scales 9. In the lower panel the 
sum of the two SFs is shown, revealing both a 
break and a shallow spectrum, as observed in Fig- 
ure 3. 

Due to a larger electron density and magnetic 
field strength in the spiral arms (Beck ct al. 1996; 
Frick et al. 2001; Cordes & Lazio 2003), the ma- 
jor contribution to the RM occurs in the spiral 
arms. Gaensler et al. (2001) determined that the 
observed polarized emission originates predomi- 
nantly from the Crux arm at a distance of 3.5 kpc 
from the Sun. The major contribution to the RM 
is most likely made by the Carina and Local spiral 
arms, located at 100 pc and 1.5 kpc from the Sun, 
respectively. Then, the Local arm corresponds to 
the nearby screen in Figure 6, and the Carina arm 
to the far screen. The break in the SF can be inter- 
preted as the location of the outer scale of struc- 
ture in the far screen, and therefore denotes the 
outer scale of the fluctuations in the Carina spiral 
arm with rfar.o ~ 2 pc. Similarly, a lower limit 
to the outer scale of fluctuations in the nearby 



screen is given by the position at which the large- 
scale anisotropy starts dominating the spectrum, 
i.e. log(r) = 0. This corresponds to the outer scale 
of structure in the Local Arm rnear.o ^ 1-7 pc. 
The shallow slope of SFem conflrms this hypothe- 
sis, as the Ha emission in this direction is believed 
not to originate farther away than 1 2 kpc, so 
all Ha emission from the Crux arm would be ab- 
sorbed by intervening dust. Note that Enfilin & 
Vogt (2003) demonstrated that the autocorrela- 
tion length of the Galactic magnetic field is smaller 
than the autocorrelation length of the RM in typ- 
ical astrophysical situations, so that the given val- 
ues for the outer scale of structure could be an 
upper limit. 

5. Hii REGIONS AS THE DOMINANT 
SOURCE OF STRUCTURE 

The only other comparable study of EM and 
RM fluctuations is that made by Minter & Span- 
gler (1996). Minter & Spangler found an outer 
scale for three-dimensional turbulence in ionized 
gas of around 4 pc in a region of about 100 square 
degrees centered at {l,b) = (143°, —21°), but con- 
cluded that two-dimensional turbulence extended 
to larger scales of at least several hundred parsecs. 
Similarly, Simonetti & Cordes (1986) found struc- 
ture in the RM of extragalactic point sources in 
a region at 70° < I < 110°, -45° < 5 < 5° on 
scales of 4° up to 40°, equivalent to 35 to 400 pc 
at a distance of 2 kpc. Studies of velocity fluc- 
tuations in neutral gas also provide evidence for 
turbulence out to scales of 100 - 1000 pc (Larson 
1979; Padoan et al. 2001). 

These results are all in striking contrast to the 
structure function analyses presented in Section 3, 
where we provide evidence that the outer scale for 
fluctuations in the density of ionized gas is only 
To = 1 — 2 pc. Because of the ubiquity of the turbu- 
lent cascade otherwise seen in ionized gas over at 
least 12 orders of magnitude in scale (Armstrong 
et al. 1995), we are reluctant to interpret our result 
as indicating that a different physical mechanism 
is operating in our particular region. Rather, we 
interpret our results as being due to an additional 
component of turbulence in the warm ISM. We 
propose that there are two contributions to tur- 
bulence in ionized gas: the overall cascade from 
the very largest Galactic scales described by Arm- 
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strong et al. (1995) found thronghont the ISM; 
and a more localized source of turbulence with an 
outer scale of ~ 2 pc, possibly of much stronger 
amplitude, found only in the Galactic plane. 

There is already good evidence to support such 
a possibility. The previous analyses of fluctuations 
in Ha emission and in Faraday rotation which have 
argued for an outer scale Vq > 100 — 1000 pc, have 
focused on sources at high Galactic latitudes, in 
the outer Galaxy, or in specific localized regions. 
However, studies of turbulence including the ion- 
ized gas in the Galactic plane of the inner Galaxy 
generally show enhanced turbulence in those re- 
gions. Clegg et al. (1992) considered RMs of ex- 
tragalactic sources in the inner Galaxy and at low 
latitude (45° < I < 93°, |6| < 5°). Because their 
background sources provided only sparse and ir- 
regular sampling, Clegg et al. were unable to cal- 
culate a continuous structure function with a well- 
defined slope, as we have done here. However, 
their data provide clear evidence that at scales 
> 1°, the amplitude of RM fluctuations at low lat- 
itudes is a factor > ID'' stronger than those seen 
at high latitudes. Clegg et al. concluded that this 
results from enhanced turbulence at low latitudes, 
generated by discrete structures in the Galactic 
plane. Simonetti & Cordes (1986) also argued for 
an additional contribution to the turbulence in the 
Galactic plane by comparing RMs of extragalactic 
sources in and out of the Galactic plane. 

A similar conclusion was reached by Spangler 
& Reynolds (1990), who compared the EMs to- 
wards eight extragalactic sources with the sizes 
of these sources' scattering disks as measured by 
VLBI. They found that heavily scattered sources 
were the same sources which also showed enhanced 
EMs along their sight lines. They concluded that 
heavily scattered sources are viewed through an 
additional component of ionized gas, with turbu- 
lent properties distinct from those of the diffuse 
ionized medium through which weakly scattered 
sources are seen. Spangler & Reynolds suggested 
that the additional component is individual Hll 
regions in the Galactic plane. Furthermore, stud- 
ies of interstellar scintillation of pulsar signals and 
angular broadening of extragalactic radio sources 
indicate higher electron density fluctuations in the 
inner Galaxy than in the solar neighborhood (Rao 
& Ananthakrishnan 1984; Dennison et al. 1984; 
Cordes, Weisberg & Boriakoff 1985). 



Several studies have tried to characterize the 
nature of this additional component of structure in 
the inner Galaxy. Cordes et al. (1985) concluded 
from pulsar scattering measurements that the ion- 
ized plasma had both clumped and nearly uniform 
components. They suggested that the enhanced 
scattering in the clumped component could be due 
to ionization fronts associated with H II regions or 
shocks associated with stellar winds and super- 
nova shells. Three different components have been 
identifled by Pynzar' & Shishov (1999), based 
on observations of pulsar scintillation timescales, 
pulse broadening and dispersion measures. They 
distinguish a component A as the diffuse gas at 
high latitudes, which is statistically uniform; a 
clumpy component BI which is distributed uni- 
formly in the spiral arms and may correspond to 
Stromgren spheres of 07 to BO stars; and a clumpy 
component BII which is concentrated in "com- 
pact regions" (~ 1.5 kpc) and is associated with 
known H II regions or supernova remnants. A sim- 
ilar subdivision was made by Ehle & Beck (1993) 
for the external spiral galaxy NGC 6946. Ehle 
& Beck created a model explaining Faraday ro- 
tation, depolarization and thermal radio emission 
as observed at various wavelengths from A2.8 cm 
to A20.5 cm. Their best fit model contains ion- 
ized gas in three components: (1) the diffuse ion- 
ized gas with low density and high filling factor; 
(2) the classical giant Hll regions with high den- 
sity and low filling factor; and (3) small (^ 1 pc) 
low-density H ii regions with intermediate filling 
factor. 

Our measurements appear to confirm this ear- 
lier work, in that we similarly have identified a re- 
gion for which the characteristics of electron den- 
sity fluctuations arc distinct from those found at 
high latitudes and in quiescent regions. While Si- 
monetti & Cordes (1986), Spangler & Reynolds 
(1990) and Clegg et al. (1992) showed that the 
amplitude of density fluctuations was enhanced 
on parsec scales in the inner Galaxy, we have 
demonstrated that the outer scale of such fluc- 
tuations is also widely diflerent. We note that 
Spangler & Reynolds suggested that the param- 
eter which differentiates the two types of turbu- 
lence is o? jro, where a is the standard deviation 
of the density fluctuations divided by the rms den- 
sity, and represents the depth of density modula- 
tion due to turbulence. Spangler & Reynolds es- 
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timatcd Tq ~ 180 pc for diffuse ionized gas, 
and To « 1.0 a'^ pc for regions of enhanced scat- 
tering. Assuming strong turbulence (a w 1) for 
the regions of strong scattering, they found a typ- 
ical scale of 1 pc, which is very similar to what 
we observe here. Other indications that an outer 
scale of stnicture of a few parsec may exist in 
the Galactic plane are presented by Haverkorn et 
al. (2003a). They show SFs of RM from the dif- 
fuse synchrotron background in two regions in the 
second quadrant at latitudes of 8° and 16°. Be- 
cause they determined the RM from observations 
at 350 MHz, only the nearby medimn 500 pc) is 
probed, i.e. only the thin disk. The SFs computed 
by Haverkorn et al. may exhibit a break at scales 
of about 2 pc, and are flat on larger scales. Fur- 
thermore, the SFrm of extragalactic sources pre- 
sented by Lazio, Spangler & Cordes (1990) seem 
to saturate on scales of a few parsec, although they 
do not mention this feature. Finally, Simonetti & 
Cordes (1986) discuss RMs from a region in the 
Galactic plane and show that there is no structure 
in this region on scales above 4°. 

We therefore confirm the supposition made by 
earlier authors that in the denser and more compli- 
cated regions of our Galaxy found at low latitudes, 
the dominant source of injection for turbulence 
appears to be relatively small, discrete, sources. 
The density fluctuations induced by these sources 
dominate the turbulent cascade generated at much 
larger scales which is seen at higher latitudes. 

What are the sources likely to be injecting such 
turbulence? Ionized gas in the Galactic plane 
is dominated by the contributions from individ- 
ual H II regions from massive stars. While such 
sources span a very wide range of sizes, depending 
on the mass and number of the central powering 
stars and on the ambient density, there is reason- 
able evidence that the characteristic scale for such 
sources is 1 2 pc, as implied here. The radius 
of an idealized Stromgren sphere for the late-type 
B stars which most likely dominate the population 
of photo-ionizing sources in the Galaxy is indeed 
a few pc (e.g. Prentice & ter Haar 1969), while 
measurements of the sizes of extragalactic H ii re- 
gions show that the distribution indeed peaks at 
diameters < 5 — 10 pc (Hodge, Lee & Kennicutt 
1989; Paladini, Davies & DeZotti 2003). 

If the observed structure is mainly due to H ii 
regions from B stars, their number should be suffi- 



cient to cover the entire field of view. We can esti- 
mate the number of B3 to AO stars per square par- 
sec of the disk from the Present Day Mass Func- 
tion (PDMF) presented by Miller & Scalo (1979), 
and the size of their Hii regions from Prentice 
& Ter Haar (1969) (who assume a mean density 
n ~ 1 cm""^). Stars more massive than B3 are ex- 
cluded because of their scarceness, and stars later 
than AO because their H ii regions are of negligible 
size. We construct a three-dimensional model of 
our field of view, placing B3 to AO stars randomly 
in the Local and Carina spiral arms according to 
the PDMF, and no stars anywhere else. We as- 
sume a thickness of 160 pc for both arms (twice 
the scale height) and the Sun positioned at the 
far side of the Local arm. The resulting distri- 
bution of Hii regions is shown in Figure 7. The 
upper plot shows a two-dimensional cut through 
the model volume, where the x-axis is directed 
along the line of sight, the observer is located at 
(x, y) = (0, 0) and the dashed lines give the field of 
view of the part of the SGPS Test Region used for 
the structure function calculation. The lower left 
plot shows the distribution of the sizes of H II re- 
gions, and the lower right plot presents how much 
of the line of sight is occupied by Hii regions, for 
all independent beams in the region. Clearly, the 
number of H II regions is enough to dominate the 
amount of structure at every line of sight. 

Interestingly, we note that Armstrong et al. 
(1995) calculated the contribution of the warm 
ionized ISM to interstellar turbulence, and con- 
cluded that spheres of density rig ~ 0.2 cm~^ and 
radius ~ 2 pc could produce fluctuations compara- 
ble to those observed. Armstrong et al. suggested 
that the amplitudes of fluctuations on such scales 
were consistent with the turbulent cascade seen at 
both larger and smaller scales. A calculation of 
whether our data contain evidence for higher am- 
plitude than or comparable fluctuations to those 
found elsewhere will be the subject of a later paper 
in this series. 

Mac Low & Klessen (2004) argue that the 
amount of energy input from H ii regions is only 
a few percent of the energy needed to drive the 
turbulence in diffuse ionized gas. However, the 
additional component of turbulence that we ob- 
served could be almost entirely due to structure 
within the H ii regions themselves, with a negligi- 
ble contribution from the turbulence in the diffuse 
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Fig. 7. — Top: two-dimensional cut through a volume distribution of H II regions of B3 to AO stars. The 

observer is located at (x, y) = (0, 0), and the dashed lines indicate the field of view of the part of the SGPS 
Test Region used for the structure function analysis. Stars are randomly positioned in the Local and the 
Carina spiral arms. The circles indicate the sizes of the individual Stromgren spheres. Bottom: number of 
Hil regions per cubic parsec as a function of their radii (left) and the amount of the line of sight taken up 
by Hii regions for many lines of sight in the field of view (right). 
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warm ISM (Spangler & Reynolds 1990; Clegg et 
al. 1992). 

If individual Hll regions dominate the struc- 
ture, what is the expected SF spectral index 5? 
The additional structure component in the ionized 
gas could be due to discrete edges of Hii regions, 
turbulence inside the Hii regions, turbulence in 
the ISM invoked by the Hll regions, or a combi- 
nation of these. For structure in RM caused only 
by discrete edges of H ii regions SFrm is expected 
to exhibit a slope 6=2 (Rickett, private commu- 
nication). There is ample observational evidence 
that Hii regions arc turbulent (O'Dell 1991; Jon- 
cas 1999; Gaensler et al. 2001), but this turbulence 
may operate on much smaller scales of ~ 0.1 pc 
(Joncas 1999). The observed steepening of the 
SFrm slope below log(r) « —1.1 possibly consti- 
tutes a transition region to steeper slopes of 6 = 2, 
or to Kolmogorov-like turbulence {b = 5/3) on 
scales < 2 pc. 

6. CONCLUSIONS 

We have analyzed structure functions (SF) of 
RM and EM data in the SGPS Test Region. These 
show very consistent results: 

• the SFs of both RM and EM, SFrm and 
SFem, exhibit a linear slope in log- log space, 
with SFrm showing a break at ~ 4' to a 
flatter slope at larger angular scales; 

• the slope at angular scales > 4' is shallow 
though non-zero in both SFrm and SFemi 
with a SF spectral index 6 w 0.2; 

• at scales larger than about a degree, the 

spectral index of SFrm and SFem is anisotropic, 
forming a quasi-sinusoidal dependence on 
position angle. 

The anisotropy in one- dimensional SFs at large 
scales is explained by a large-scale gradient in elec- 
tron density, possibly accompanied by magnetic 
field, across the field of view due to a foreground 
structure. The break in SFrm and shallowness 
of the slope in both SFrm and SFem at smaller 
scales can be explained by the superposition of 
two contributions to EM and RM of the Local 
and Carina spiral arms. Within such a model, 
we infer the outer scale of structure in the spi- 
ral arms to be about 2 pc. As Kolmogorov-like 



turbulence is observed in the ISM on scales much 
larger than a few pc, these results constitute evi- 
dence for an additional contribution to turbulent 
fluctuations in the Galactic plane. The inferred 
outer scale of fluctuations agrees with the size 
of Stromgren spheres around the most abundant 
late-type B stars. Therefore, Hii regions could 
provide the dominant source of structure on pc 
scales in the Galactic spiral arms. 

The suggestion that the break in SFrm and the 
shallowness of both SFrm and SFem is due to two 
separate Faraday screens will be tested in a forth- 
coming paper by studying SFs in different regions 
of the sky. This can also shed light on the oc- 
currence of any spatial variations of the turbulent 
spectrum. Furthermore, the combined RM and 
EM data should enable at least partial decoupling 
of magnetic field and electron density, and the am- 
plitudes of the SF will yield information on this 
additional component of structure in the ionized 
ISM in the Galactic plane. 
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A. NOISE CONTRIBUTION TO THE STRUCTURE FUNCTION 

For a rotation measure RM ± 6 rad m"^, the SFs of the RM and the error d can be separated using 

SF(RMo6,) = ([(RM(a;) + S{x)) - {RM{x + r)+ S{x + r))f) 

= (RW^ix) + RW^ix + r) - 2RM(a;)RM(a; + r)) + 

(^^(a;) + 6'^{x + r)- 25{x)5{x + r)) 
= SF(RM) + 2 (,^2^ (Al) 
= SF(RM) + SF((5) (A2) 

since (RM(5) = and {5{x)5{x + r)) = 0. Therefore, the contribution of noise in RM (which is assumed 
Gaussian) can be taken into account by subtracting a SF of a Gaussian distribution of noise with width 
from the SF of the observed RM. 
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